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Abstract 



We have searched for a neutral H dibaryon decaying via H — > An and 
H — > Our search has yielded two candidate events from which we set an 
upper limit on the H production cross section. Normalizing to the inclusive A 
production cross section, we find {da H /ain) / (da A /dn) < 6.3 x 1(T 6 at 
C.L., for an H of mass «2.15 GeV/c 2 . 
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The theory of quantum chromodynamics imposes no specific limitation on the number 
of quarks composing hadrons other than that they form color singlet states. Although only 
qqq and qq states have been observed, other combinations can form color singlets. Jaffe |]J 
has proposed that a six-quark state uuddss may have sufficient color-magnetic binding to be 
stable against strong decay. Such a state, which Jaffe named H, would decay weakly, and 
the resultant long lifetime would allow the possibility of observing such particles in neutral 
beams. Theoretical estimates of m H have varied widely, ranging from a deeply bound 
state with m H < 2.10 GeV/c 2 to a slightly unbound state with m H near the AA threshold, 
2.23 GeV/c 2 . In this mass range the H would decay almost exclusively to An, S ra, and 
||. Several previous experiments have searched for if's but with no compelling success 
Q. The search described here is sensitive to if's having mass and lifetime in a previously 
unexplored range. 

We have searched for H —>■ An and H — > S°n — > A'jn decays by looking in a neutral beam 
for A — > pn~ decays in which the A momentum vector does not point back to the production 
target. The experiment, E888, was performed in the B5 beamline of the Alternating Gradient 
Synchotron (AGS) of Brookhaven National Laboratory. A second phase of the experiment 
searched for long-lived fTs by using a diffractive dissociation technique 0. The detector 
used for the decay search (Fig. |IJ) was essentially that used for the E791 rare kaon decay 
experiment and has been described in detail elsewhere 0. In brief, a neutral beam was 
produced using the 24 GeV/c proton beam from the AGS incident on a 1.4 interaction length 
Cu target. The targeting angle was 48 mr. After passing through a series of collimators and 
two successive sweeping magnets, the neutral beam entered a 10 m long vacuum decay tank 
within which candidate A's decayed. Downstream of the tank was a two arm spectrometer 
consisting of two magnets with approximately equal and opposite p T impulses and 5 drift 
chamber (DC) stations located before, after, and in between the magnets. Downstream of 
the spectrometer on each side of the beam were a pair of trigger scintillator hodoscopes 
(TSCs), a threshold Cherenkov counter (CER), a lead-glass array (PbG), 0.91 m of iron 
to filter out hadrons, a muon-detecting hodoscope (MHO), and a muon rangefinder (MRG) 
consisting of marble and aluminum slabs interspersed with streamer tubes. For the first 
half of the run the Cherenkov counters were filled with a He-N mixture (n = 1.000114) to 
identify electrons; for the second half the left-side counter was filled with freon (n = 1.0011) 
to identify protons from A — > pn~ (due to lack of light). Only the left counter was used for 
this purpose as the soft pion from A — > pir~ decay is accepted only when on the right; when it 
is on the left, the first magnet bends it back across the beamline and it is not reconstructed. 
The lead-glass array (PbG) consisted of two layers: a layer of front blocks 3.3 radiation 
lengths (r.l.) deep and a layer of back blocks 10.5 r.l. deep. The PbG was used to identify 
electrons by comparing the total energy deposited (E tot ) with the track's momentum. A 
minimum bias trigger was defined as a coincidence between all 4 TSC counters and signals 
from the 3 most upstream DC stations. A Level 1 trigger (LI) was formed by putting 
minimum bias triggers in coincidence with veto signals from the Cherenkov counters and 
muon hodoscope. All events passing LI were passed to a Level 3 software trigger which used 
hit information from the first 3 DC stations to calculate an approximate two-body mass. 
Events with m p7r _ < 1.131 GeV/c 2 were written to tape. 

Offline, all events containing two opposite-sign tracks forming a loose vertex were kine- 
matically fit || and subjected to the following cuts: there could be at most one extra 
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track-associated hit or one missing hit in the ten DC planes which measure the x (bending) 
view of each track; the x 2 's per degree of freedom resulting from the track and vertex fits had 
to be of good quality; the A vertex had to be within the decay tank and downstream of the 
fringe field of the last sweeper magnet; both tracks had to be accepted by CER, PbG, MHO, 
and MRG detectors and have p> 1 GeV/c; neither track could intersect significant material 
such as the flange of the vacuum window; to reject background from K\ — > 7r°7r + 7r~, m n+n - 
had to be > (tti Kl ~ m n°)'> an< ^ to re j ec t background from K s — > 7r + 7r~ resulting from sec- 
ondary interactions, \m^ +n - —m KL \ had to be > 4 times the mass resolution of K\ — > 7r + 7r~ 
decays (1.55 MeV/c 2 ). 

Events passing these cuts were subjected to particle identification criteria in order to 
reject background from K\ — > -nev (K e3 ) and K\ — > nfiu (K 3 ) decays. To reject electrons, 
we require that there be no track-associated Cherenkov hit and that tracks with p > 2 GeV/c 
(< 2 GeV/c) have E tot /p < 0.60 (< 0.52). The low momentum track on the right side of 
the detector was required to deposit < 0.66 x E tot in the front PbG blocks. To reject muons 
which passed the MHO veto in the trigger, we cut events with a hit in the MRG which 
was consistent with the projection of a track and which corresponded to at least 65% of the 
expected range of a muon with that track's momentum. 

Lambda candidates were selected by requiring that \m p7T - — m A \ be less than 4 times 
the mass resolution of A — > pir~ decays (0.55 MeV/c 2 ). The data was then divided into 
two streams: a normalization stream consisting of A's which project back to the production 
target, and a signal stream consisting of A's which do not. The former were selected by 
requiring that the square of the collinearity angle # A be less than 1.5 mrad 2 , where # A is the 
angle between the reconstructed A momentum vector and a line connecting the production 
target with the decay vertex. This sample contains negligible background. The signal sample 
was selected by requiring that p T > 145 MeV/c, where p T is the A momentum transverse to 
the line connecting the production target with the decay vertex. This cut value was chosen 
to eliminate S° — > Air decays, which have a kinematic endpoint of 135 MeV/c. The p T 
distribution of A's from two-body H — > An decays exhibit an approximate Jacobian peak 
(not exact because the vertex is the A's) with an endpoint which depends upon m H . A large 
fraction of high-p T A's were found to project back to a collimator located just upstream of 
the decay tank. We thus required that the point in our beamline to which a A projects back 
be located downstream of this collimator: z pro j > 9.65 m. 

A signal region for H candidates was defined by the criteria p T > 174 MeV/c and 
N T > 5, where N T is the distance in proper lifetimes between the decay vertex and the 
nearest material (beamline element) to which the momentum vector projects back. The 
p T cut rejects K e3 decays which survive the CER, PbG, MHO, and MRG vetoes due to 
detector inefficiency, while the N T cut rejects A's which originate from collimators, flanges, 
and other beamline elements. All cuts were determined without looking at events in the 
signal region, in order that our final limit on if's be unbiased. After fixing cuts we looked 
in the signal region and observed two events. The estimated background is 0.15 events 
from A's originating from beamline elements, and < 0.21 events from K t3 decays (all K e3 
as the p T is too high for K 3 ). The former is estimated by studying the N T distribution of 
A's originating from a "hot" flange located immediately upstream of 9.65 m. The latter is 
estimated by first counting the number of final events cut because the low-momentum track 
had E tot /p > 0.7 (these are electrons); this is then multiplied by the ratio of the number 
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of electrons passing PbG analysis cuts to the number having E tot /p > 0.7, as determined 
from a sample of K e3 decays. The N T vs. p T plot for the final high-p T A sample is shown in 
Fig. |2|. In this figure the Cherenkov veto for the freon counter is not imposed. A band of 
K^ 3 decays is visible at p T ~ 150 MeV/c which results from the p T > 145 MeV/c cut and 
the m wn > (m K — m^o) cut; this latter cut constrains p T from above. A's which originate 
from beamline elements are visible at low N T . For the freon subset, when we require that 
there be no signal in the left Cherenkov counter, all but two K 3 decays are eliminated while 
all A's at low N T remain. 

Also visible in Fig. ^| are our two candidates, which have p T of 187 and 191 MeV/c and 
N T of 6.7 and 9.4. The p T values correspond to a Jacobian peak from H — ► An decay if 
m H pa 2.09 GeV/c 2 . The probability for a K 3 decay to have such high p T is extremely 
small, as it is kinematically forbidden for a K 3 decay to have both > [jn Kh — m^o) and 
p T > 160 MeV/c (Fig. |3|). The probability for a K e3 decay to look like these events is also 
very small, as the PbG response for the electron candidate tracks is very uncharacteristic 
of electrons: E to Jp = 0.44 and 0.27, and for both events E {ront /E tot = (Fig. |). This 
response is typical of pions from A — > pn~ decay. To investigate background from neutrons 
in the beam interacting with residual gas molecules in the decay tank, we recorded and 
analyzed a sample of data equivalent to 1% of the total sample with the decay tank vacuum 
spoiled by a factor 2.7 x 10 3 . This sample yielded one event in the signal region, implying a 
background level in the rest of the data of 0.04 events. We also studied potential background 
from S° — > A7T° decays where the S° originates from a beamline element; from Monte Carlo 
simulation and the number of A's observed originating from beamline elements, we estimate 
a background of less than 0.10 events. The total background estimate from known sources is 
less than 0.50 events. The probability of 0.50 events fluctuating up to two or more events is 
0.090; if such a fluctuation occurred, it is remarkable that the p T of the events is so similar. 

A 90% C.L. upper limit on the H production cross section can be expressed in terms of 
the inclusive A production cross section as follows: 

da H j A A B(A -> pii~) da A 

dn < N A rg A H B(H -> AX) dQ ' U 

where N A rg = 20 433 is the number of A's originating from the target, A A and A H are 
geometric acceptances for A's originating from the target and from H decays, respectively, 
B(A — > pn~) and B(H — > AX) are branching ratios, da A /dQ is the inclusive A production 
cross section, and £ is the factor which multiplies the single-event sensitivity to give the value 
of dan/dQ which has a 10% chance of producing < 2 detected events. Here we conservatively 
assume no background and take £ = 5.32. The acceptance A H accounts for the fact that A's 
from H's must project back to a restricted region of the beamline. Since A — > pn~ decays 
are common to both signal and normalization channels, all trigger and detection efficiencies 
divide out of Eq. (P . 

The acceptances A A and Ah were determined from Monte Carlo simulation using several 
different estimates of the production momentum spectra. For the H simulation, a central 
production spectrum was used with a broad peak at x F = 0. A spectrum corresponding 
to a AA coalescence model for H production |7| resulted in a limit on dan/dVL about 50% 
lower. We quote here the more conservative limit resulting from the central production 
spectrum. The inclusive A production spectrum was taken from a measurement by Abe et 
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al. ||; comparison with our data shows very good agreement. 

The acceptance Ah also depends crucially on H lifetime. Here we assume the relationship 
between m H and t h , B(H — > S° n), and — > A n) calculated by Donoghue et al. 0, and 
obtain 90% C.L. upper limits on (dan/dQ) / {da^/dVl) as a function of m^. Our acceptance 
is maximum for r H « 8 ns and becomes small for r fl ^ 1 ns due to the z j > 9.65 m cut. 
Our limits for (dan/dfl) / (da\/dQ) are plotted in Fig. ||. For m H « 2.15 GeV/c 2 , Jaffe's 
original prediction, 



da 



H 



dn 



48 mr 



< 6.3 x 10" 



do* 



dn 



(90% C.L. 



48 mr 



From Abe et al. 



da\/dQ\ 



48 mr 



366 mb/sr, so dan/dfl\ 



48 mr 



(2) 



< 2.3 /ib/sr. For 



= 2.09 GeV/c , consistent with the observed A p T , the acceptance is lower and the 
two candidate events correspond to a differential cross section of 44 ^ /xb/sr. The authors 
of Ref. |3| note that t h may be shorter than their predicted value by up to a factor of two; 
this would increase our acceptance for m H ^ 2.18 GeV/c 2 and decrease our acceptance 
for m H greater than this value. The resultant 90% C.L. upper limits are plotted as the 
dashed contour in Fig. |5]. If we assume that the invariant cross section Ed 3 a/dp 3 has the 
form A(l — \x\) e p t : then our limit (H) corresponds to o u < 60 nb for a wide range of 
parameters B and C. 

There are few theoretical predictions of the H production cross section. Cousins and 
Klein predict a differential cross section of ^ 100 /xb/sr for p-Cu interactions at our 
targeting angle based on a AA coalescence model. Cole et al. considers AA and S°n 
coalescence and predicts cr tot ~ (3 x 10~ 5 ) x cr inelas for p-Cu collisions at AGS energies; 
taking the inelastic cross section a inelas to be ~ 780 mb [0] gives cr tot ~ 23 /ib. Rotondo 



1 L| considers only 2° n coalescence and predicts a total cross section at Fermilab energies 



of 1.2 /ib. 
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FIG. 1. The E888 detector and beamline. 
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FIG. 2. N T vs. p T for the high-p T A sample. The signal region is denoted by dashed lines. The 
band of events from p T = 145 to p T w 150 MeV/c are K^ 3 decays; the leftmost edge is due to a p T 
cut, while the rightmost edge is due to a lower cut on m n7T . 
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FIG. 3. m wn vs. p T for the final high-p T A sample. The two events in the signal region are 
circled. The cluster of events at p T ~ 150 GeV/c, m n7T ~ 365 GeV/c 2 are consistent with Monte 
Carlo simulated K 3 decays. There is a third event which is well-separated from the K 3 decays 
and which lies just outside the signal region; the 3 separated events are consistent with H — ► An 
decay if m H « 2.09 GeV/c 2 . 
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FIG. 4. Ef ront /E t ot (PbG) vs. E tot /p for: a) the low momentum track of A's from the final 
high-p T sample, and b) low momentum electrons from K e3 decay. In (a), the tracks from the two 
events in the signal region are circled. There are 4.7 times as many events in (b) as in (a). 
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FIG. 5. 90% C.L. upper limits on the H production cross section as a function of m H or t h 
(see Ref. [3]). The dashed contour corresponds to an H lifetime half that given on the top scale. 
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